arXiv:1505.02876vl [hep-lat] 12 May 2015 


Transition of p —)■ 7r7 in Lattice QCD 


Benjamin J. Owen/ Waseem Kamleh/ Derek B. Leinweber/ M. Selim Mahbub/’^ and Benjamin J. Menadue^’^ 

^Special Research Centre for the Subatomic Structure of Matter (CSSM), 

Department of Physics, University of Adelaide, South Australia 5005, Australi?^ 

^Digital Productivity Flagship, 

CSIRO, 15 College Road, Sandy Bay, TAS 7005, Australia 
^National Computational Infrastructure (NCI), 

Australian National University, Australian Capital Territory 0200, Australia 

(Dated: May 13, 2015) 

With the ongoing experimental interest in exploring the excited hadron spectrum, evaluations 
of the matrix elements describing the formation and decay of such states via radiative processes 
provide us with an important connection between theory and experiment. In particular, determina¬ 
tions obtained via the lattice allow for a direct comparison of QCD-expectation with experimental 
observation. Here we present the first light quark determination of the p —>■ TT'y transition form 
factor from lattice QCD using dynamical quarks. Using the PACS-CS 2-1-1 flavour QCD ensembles 
we are able to obtain results across a range of masses, to the near physical value of m-,, = 157 MeV. 

An important aspect of our approach is the use of variational methods to isolate the desired QCD 
eigenstate. For low-lying states, such techniques facilitate the removal of excited state contributions. 

In principle the method enables one to consider arbitrary eigenstates. We find our results are in 
accord with the non-relativistic quark model for heavy masses. In moving towards the light-quark 
regime we observe an interesting quark mass dependence, contrary to the quark model expectation. 

Comparison of our light-quark result with experimental determinations highlights a significant dis¬ 
crepancy suggesting that disconnected sea-quark loop contributions may play a significant role in 
fully describing this process. 


I. INTRODUCTION 

In the search for a complete understanding of the 
hadron spectrum, hadronic transitions provide a crucial 
dialogue between experiment and lattice QCD. Experi¬ 
mental programs such as those at JLAB and Mainz are 
producing large bodies of high-quality data [I| provid¬ 
ing valuable insight into the underlying structure and 
dynamics of hadronic excitations. Meanwhile, signifi¬ 
cant effort by the lattice community has seen remarkable 
progress in understanding and mapping out the hadron 
spectrum from the perspective of QCD I, i. Hav¬ 
ing isolated such states, evaluations of their transition 
form factors provide crucial insight into their produc¬ 
tion mechanisms allowing for direct comparison between 
QCD and the experimental observations. Another excit¬ 
ing prospect is that lattice determinations may help in 
the experimental search for as yet unseen QCD eigen¬ 
states, such as exotic hadron states being sought by the 
GlueX experiment Q. 

Among the simplest radiative process that can be ex¬ 
amined on the lattice is the radiative decay of the rho me¬ 
son into a pion. This process is described by a single form 
factor, Gmi{Q^), and as both states are the lowest lying 
single particle states in their respective channels, extrac¬ 
tion is relatively simple. Evaluations of this transition 
were considered early in the course of lattice structure 
calculations, first by Woloshyn Q and later by Lubicz 
and Crisafulli More recently, our preliminary study 
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0 examined both the p and K* transition form factors 
down to lighter masses (m^ ~ 300 MeV). However, all of 
these studies used quenched gauge field configurations. 
Currently the only evaluations in full QCD are those of 
Edwards @ and Shultz Q in which this transition along 
with a range of other elastic and transitions form fac¬ 
tors were considered. However, the parameters consid¬ 
ered were far from physical, particularly the large pion 
mass. There is thus a need to examine this transition at 
light quark masses using dynamical quarks. 

Recently the issue of excited states has received re¬ 
newed focus dMl as lattice determinations of several 
quantities, most notably the nucleon axial charge, gA, 
and quark momentum fraction, (x), continue to differ 
systematically from the experimental results despite op¬ 
erating at near physical quark masses and large lattice 
volumes. In Refs, [l^, Hjli presented a framework 
for utilising variational techniques for the extraction of 
hadronic matrix elements. For both the calculation of gA 
0 and the electromagnetic form factors of the tt and p 
mesons M , we found that the use of correlation ma¬ 
trix methods led to significantly improved ground state 
dominance over the standard single smeared source ap¬ 
proach. This in turn allowed for the use of earlier current 
insertion times and subsequently consider earlier fit win¬ 
dows resulting in improved statistical uncertainties. Fur¬ 
thermore, within Ref. 0 we were able to examine the 
corresponding matrix elements for both excited pi and 
rho meson states. The same techniques were utilised by 
Shultz et al. Q to examine the transition form factors 
for the lightest tt and p excited states. In this work we 
will make use of such techniques in our determination of 
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the p —>■ 7r7 transition form factors across a wide range 
of quark masses including the light-quark regime. 

The remainder of this investigation is organized as fol¬ 
lows. In section 2 we present a brief summary of the 
use of variational methods in the evaluation of transition 
matrix elements. Section 3 outlines how we intend to ex¬ 
tract the transition form factor, G'mi(Q^) and finally a 
summary of calculation details in Section 4. Our results 
are discussed in Section 5 followed by concluding remarks 
in Section 6. 


II. CORRELATION MATRIX METHODS FOR 
HADRON TRANSITIONS 

The use of variational techniques for the calculation 
of hadronic matrix elements has been presented in Refs. 
[T^ fl^ and here we present a brief outline to define 
our procedure and notation. Similar prescriptions have 
been explored in Refs, [isl - l^ . The variational method 
is an approach to generating a set of ideal opera¬ 
tors, {0“}, constructed to have overlap with an individ¬ 
ual eigenstate 

■ ( 1 ) 

This is achieved by choosing an existing basis of opera¬ 
tors, {xi}, and constructing these idealised operators as 
linear combinations 

=Y^v^X^ix), =Y^Xjix)u'^ . ( 2 ) 

* j 

Starting from the matrix of two-point correlation func¬ 
tions 

G^M 0 = E ^ I I ^’ 

X 

one can show that both vf Gtj (p, t -\- St) and Gij (p, t -I- 
St) provide a recurrence relation from which the neces¬ 
sary vectors, vf and are the solutions of the following 
generalised eigenvalue equations 

vf Gij (p, to + St) = Vi Gij (p, to) (3a) 

Gij (p, to + St) u° = Gij (p, to) . (3b) 

With the optimised operator for the state \a,p) suitably 
defined, one can construct the corresponding correlator 
by projecting the matrix of correlators 

G{p,t;a) = {v{i^)f G,j{pG) {u{p))j ■ 

Here we note that the correlators that feed into Eq. (I3a1) 
and (Hi) are functions of the 3-momentum, p. There¬ 
fore it is necessary to evaluate eigenvectors for each mo¬ 
mentum considered. This is of particular relevance to 
three-point correlators where the incoming and outgoing 
momenta in general differ. 


To obtain the three-point correlator necessary for de¬ 
scribing a particular transition, it is a simple matter 
of projecting the relevant states with the corresponding 
eigenvectors for source and sink. In order to do this, 
one first evaluates the matrix of three-point correlation 
functions required for the transition in question 

X2,Xi 

X {n\xbAx2)oix,)xl,io)\n), (4) 

where O is the current operator through which the tran¬ 
sition takes place and the labels a and b highlight that 
the source and sink operators can correspond to different 
In this situation, one must construct ideal opera¬ 
tors for each set of quantum numbers. By performing the 
variational analysis with the bases { Xa,i } and { Xb,i }, we 
obtain two sets of eigenvectors: (fa)“, (mq)^ and {vb)f, 

{ub)j which form optimised operators of type a and b 
respectively. It is then a simple matter of projecting 
the necessary eigenvectors onto the matrix of three-point 
functions, where care is taken to ensure that the projec¬ 
tion is done with the correct momenta and operator for 
source and sink 

G‘^^{p',p,t2,ti;Aa) = 

Mp'))^ {G^o^%ip\p,t2,h) (uAm ■ 

Having acquired the projected two and three-point cor¬ 
relation functions, determination of matrix elements then 
follows in the standard way through the construction of 
a suitable ratio. For transitions, we choose to work with 
a modified form of the ratio presented in Ref. [2^ 

i?(p',p;/3,a) = 

{Go^^{p'A,t2,ti;l3,a)) (G^~^°(p,p', ^2, ti; a,/?)) 

(Gaip,t2]a)) {Gbip',t2-,P)) 

( 5 ) 

Though this requires the evaluation of both +q and —q 
sequential source technique (SST) propagators 1^. 
if one wishes to use both {I/j and {U*} configurations 
as was done in Refs. [l3, and will be done here, then 
these are required regardless of the choice of ratio . 


HI. p^Tv-y TRANSITION FORM FACTORS 

Having outlined the method for obtaining the pro¬ 
jected correlators and the corresponding ratios, we now 
consider how to isolate the form factor describing the ra¬ 
diative decay of the rho meson. In order to do so we 
make the identification that a = n, b = p and O = 
the electromagnetic current. For this discussion we shall 
work with the standard Minkowski metric. Starting from 
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the projected three-point correlation function 

X2,Xi 

( 6 ) 


we begin by inserting completeness identities between our 
operators in order to obtain the desired matrix element, 
as well as operator overlap factors. This matrix element, 
which describes the transition, can be parametrised by a 
single form factor Gmi 


I ’^'"(0) I P/3 (p', s')) = 




2^E^MEp,{p') 

GMi{Q^)e>^^^^p'sP.4ip',s'), (7) 


where we note that our choice of normalisation is that of 
Ref. 0. Furthermore we note the pion overlap can be 
expressed as 

{I I np{p)) = 

^J2 (p) 


and similarly the rho meson overlap 


{n\rpm\pp{p.s)) 


Ja/3 


V^EpJp) 


^u{p,s)z^{p) 


where e(p, s) is a spin polarisation vector which satisfies 


J2^t^^Es)€l{p,s) 



Substituting these expressions into Eq. and making 
use of the spin-sum identity we arrive at 


e--E-c if) (t2-ti)g-.Epp (p')ii 

4 (pO Ep, (pO 


Z“(p)Zf(p') 


-\-ie 




GMl{Q^)e^^^^p'sPa 


Taking the Hermitian conjugate of Eq. 0, we can obtain 
the time-reversed matrix element necessary for Gt^-j^p- 
Following the same procedure we arrive at the corre¬ 
sponding expression 


{G^^^p));:{p',p,t2,ti;P,a) = 
4 Ep (pO 


zP(p')Z:^[p) 


—le 


^p0 


GMi{Q^)e^^^^p'sPa [gr.- 


PtPv 


pp 


In order to cancel out the time dependence and the over¬ 
lap factors, Z, we require the projected two-point func¬ 
tions 

G.(p, t 2 ; a) = E I C’^> 2 ) ( 0 ) I b!) 

(Gp)..(p", t2; /?) = E ^ I ^pf (^ 2 ) <f ^ (0) I ), 


where we note that repeated indices are not summed 
over. Again applying completeness and replacing spin- 
sums, these expressions reduce to 


e--E,c(R)*2 

g^{pM) = Kip)KH p) 


‘^EKp) 


iGp)Kp'M) = - 


g-Epp (p')‘2 


^Eppip') ^ 


-ZPip')ZPHp') 


^ / P^'P 

\ 


P0 


Finally, substituting in the explicit form of both two and 
three-point correlation functions into Eq. 0, we obtain 
the expression which directly relates our ratio to the form 
factor Gmi(Q^) 


Ruip',p;l3,a) 


_ Pa Gmi(Q^)_ 

2 JE^^ (f) Epp (p') ((p'^)2 - 


( 8 ) 

For our calculation we use SST-propagators evaluated 
by fixing the current. In particular, we fix the current 3- 
momentum to g = = ±^, the current polarization 

p = 3 and the vector meson polarization to be iz = 2. 
Using these kinematics, the term e^^'^^pgPa reduces to 
(PoPi —pop'i)- We choose to refine this further by taking 
one of the particles to be at rest. In doing so, we still 
have the freedom to choose which of the particles is at 
rest giving rise to the distinct kinematics choices: p' = 
p = 0 (g = -l-'C) and p' = 0, p = ^ (g = —^), each of which 
gives distinct values of . Applying these kinematics to 
Eq. ([8|) we arrive at the final expressions used in our 
evaluation of the p —>■ Try transition form factor 


Gmi(Q') 

Gmi(Q') 


2wpg Epp{^ 

e 1^ Y 




2 nT-pg / Et^^ {q) 

e 1^ Y 


Rl{0,i-(3,a). 


IV. SIMULATION DETAILS 

For this calculation we use the PACS-CS (2+U-flavour 
dynamical-QCD gauge field configurations made 
available through the ILDG [ 2 ^. These ensembles use a 
non-perturbatively G(a)-improved Wilson fermion action 
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and Iwasaki gauge action on a 32^ x 64 lattice, with pe¬ 
riodic boundary conditions. The value /3 = 1.9 provides 
a lattice spacing of a = 0.0907 fm, yielding a physical 
box length of 2.9 fm. We have access to a total of 5 light 
quark masses, with the strange quark mass held fixed. 
The resulting pion masses range from 702 MeV down to 
156 MeV. We note that when quoting results away from 
the physical point we set the scale using the Sommer pa¬ 
rameter with To = 0.492 fm [2^. Due to the limited 
number of conhgurations at the lighter masses, we make 
use of multiple quark sources on each configuration. We 
also evaluate the correlators on the {U*} conhgurations 
through the use of the 17*-trick outlined in Refs. [13,mi. 

For the evaluation of the three-point correlators we 
follow the approach where the SST propagators are eval¬ 
uated with the current held hxed. For this we use a con¬ 
served vector current Is^l with polarisation fi = 3 and 
3-momentum transfer q = This current is inserted 

at time tc = 21 relative to the quark source at tsrc = 16. 
We note that for our choice of ratio we require both +q 
and —q, however this is also required for the ?7*-trick 
when using SST-propagators. Our error analysis is per¬ 
formed using a second-order jackknife, with the x^/dof 
for our hts is obtained through the covariance matrix. 

As was done in Ref. 0, the anti-quark contribution 
is evaluated through considerations of charge-conjugation 
in order to avoid the need to evaluate backwards propa¬ 
gating SST-propagators. In doing so one finds that the 
quark sectors for the forward and backwards propagat¬ 
ing quark fields are of equal magnitude, but differ in sign 
when exact isospin is present. 

For the variational analysis, we use a basis of local me¬ 
son operators of varying widths [s^. This is achieved 
by applying different levels of gauge-invariant Gaussian 
smearing to the quark sources and sinks. By using a va¬ 
riety of widths, the resulting optimised interpolators ap- 

S r to form nodal structures in the radial wave-function 
Thus we expect our operators to have greatest 
overlap with states that contain dominant s-wave compo¬ 
nents. As these states are strongly peaked at the origin, 
we would expect these states to have the largest over¬ 
lap with standard local operators and so the dominant 
source of excited state contamination for the standard 
single source approach. The local operators we choose to 
use are 

X^{x) = d{x)xbu{x) 

Xp,i{x) = d{x)'yiu{x), 

with four different smearing widths, allowing for the con¬ 
struction of a 4 X 4 correlation matrix. The correlation- 
matrix analysis is performed using to = 17 with St = 3 
for the three heavier masses and St = 2 for the remaining 
two lightest masses. 


1.0 
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(D 
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FIG. 1. Masses of the tt (blue) and p (red) mesons. The 
lighter data points at the far left are the corresponding PDG 
averages [I[. 


V. RESULTS 

Before we examine the transition form factor, we ex¬ 
amine the isolation of the eigenstates in question. In 
Fig. □ we present the masses of the tt and p mesons ex¬ 
tracted from our variational analysis. For the lightest 
quark mass, which is only slightly heavier than physical, 
we find that both states agree well with experimental 
measurements. Nonetheless, for the rho meson we must 
take care to ensure that the state we have isolated is in 
fact the single particle state and not a scattering state. 
In Ref. [1^ it was found that despite coupling poorly to 
local operators, scattering states can still populate the 
projected correlators. At early times, the weak operator 
coupling acts to suppress their contribution relative to 
the dominant single particle state. However if the sin¬ 
gle particle state is above threshold, at large Euclidean 
times its contribution is suppressed through the exponen¬ 
tial behaviour of the correlator rendering the scattering 
state the dominant term. There it was noted that in such 
a scenario, one could identify a second region of linear 
behaviour in the logarithm of the projected correlator. 
We shall use this signature to identify how far we can 
sample the correlator and still guarantee single particle 
dominance. 

On these ensembles the p meson at rest remains lighter 
than the nearest TTTr-scattering state, however this is not 
the case for non-zero momentum. In Ref. [T^ we per¬ 
formed a detailed check of our extracted eigenstates to 
ensure that we had in fact isolated the single particle 
state, by comparing the extracted energies with the dis¬ 
persion expectation using the mass obtained with the 
rho-meson at rest. This is also compared with the cor¬ 
responding bare TTTr-scattering states in order to identify 
where we should expect to see increased mixing between 
single and multi-particle states. In all cases we found that 
the extracted energies agreed well with the dispersion 
expectation, particularly for the kinematic arrangement 
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relevant to this work. Furthermore the decay thresholds 
were found to lie well away from the extracted energies. 

For the particular kinematics and parameters consid¬ 
ered in this work, we find that the determinations of 
Gmi with the pion at rest all lie below = 0, while 
those with the rho meson at rest lie above. This gives 
us values for Gmi on either side of = 0. In order to 
compare with experiment and quark model expectations 
we require a determination of G'mi(O). To do this we 
choose to interpolate between our extracted values using 
a monopole ansatz 

Gmi{Q‘^) = ^^2 _|_ Gmi(0). (9) 

This choice is motivated by vector meson dominance 
(VMD) arguments which suggest the form factor should 
exhibit such a behaviour in the region of low Q^. In these 
models, the pole mass is identihed as the vector meson, 
i.e. A ~ vrip. In Ref. Q, it was found that the VMD 
hypothesis faired poorly with the data, however we note 
that this study was conducted with rather large values 
of stemming from their small lattice volume. Con¬ 
trary to this, the results of Edwards Q which examined 
the transition over a range of between 0.02-0.6 GeV^, 
display behaviour consistent with this expectation. 

During the extraction of the quark sector contributions 
to the form factor, we compared the time-series for the 
ratio using the correlation matrix approach and the stan¬ 
dard single smeared source and sink correlator. For all 
masses and kinematics, we again find that the correla¬ 
tion matrix method improves the quality of the plateau 
over modest levels of smearing, however not to the extent 
observed in previous works [l3 . 0. In particular, the ra¬ 
tio sampling in the time-like region requires significantly 
more Euclidean time evolution than the corresponding 
ratio sampling the space-like region. We note that in 
this case the rho meson carries the momentum and the 
pion is at rest. Eigure [2] highlights this comparison for a 
single quark mass. 

In Eig. Owe present our results for the extracted val¬ 
ues for the w-quark sector contribution to the form factor, 
as well as the corresponding interpolations used to 
extract G\ji{Q). Here we choose to label these as the 
quark sector contributions to the positive-charge eigen¬ 
state of the corresponding iso-triplet. That is, the quark 
contribution is labelled as the w-quark sector while the 
anti-quark contribution is labelled as the d-quark sector. 
As was mentioned in the previous section, the anti-quark 
contribution is equal in magnitude with opposite sign, 
G'lji = —G\pi, and so we choose to show the quark con¬ 
tribution only. The quark sector contributions are for 
quarks of unit-charge. We note that the spread of 
sampled is much larger in the time-like region due to the 
increasing Qo-component stemming from the Goldstone 
nature of the pion. In the space-like region, we see this 
effect is suppressed when the momentum is coupled to 
the pion. This tight grouping of values shows a clear 
decrease in the value of as the quark masses get 



(a) with the tt at rest. 



(b) with the p at rest. 

FIG. 2. A comparison of the u-quark contribution to the 
transition form factor, G'hi, as a function of Euclidean sink 
time for a single level of smearing and our variational ap¬ 
proach. The two data sets are offset for clarity. The upper 
figure is for the tt meson at rest with the (blue) circles de¬ 
noting the results from the variational approach while the 
(red) squares illustrate traditional results using the standard 
single-source method with a moderate level of smearing (35 
sweeps). The lower figure is for the p meson at rest with 
the (green) triangles denoting the variational method and the 
(purple) diamonds denoting the single-source method. The 
vertical dashed line indicates the position of the current in¬ 
sertion. The fitted value from the variational approach has 
been included (shaded band) to highlight where the single 
source approach is consistent with our improved method. 


lighter. 

In the time-like region, the most striking feature is the 
significantly small value obtained for the lightest mass. 
However, examination of the projected two-point func¬ 
tion appears consistent with single particle state. We 
shall note that the plateau for this particular extraction 
of G'^t^ differed in nature to those at heavier masses. 
Based upon the one is able to fit much earlier, 
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FIG. 3. Quark sector results for the transition form factor, 
Gjvfi- Each curve corresponds to a different value of m,n-, 
with the top curve (blue) corresponding to the heaviest and 
masses getting lighter as we move down. For each mass we 
have access to two values of stemming from the freedom 
to choose which state is as rest. The solid line and coloured 
bands are the resulting monopole parametrisation extracted 
from the two data points, allowing us to access Gmi(0). We 
also include a VMD estimate, dashed line, obtained from us¬ 
ing the right positive data only, with p-meson mass used 
as the monopole mass A. 



FIG. 4. Monopole mass. A, (blue) obtained from applying 
monopole ansatz to the data. The corresponding p-meson 
mass (red) is included to test the VMD hypothesis. The 
monopole mass for the lightest quark mass is not included 
due to the large uncertainty in A^ stemming from the similar 
values for Gmi (QI) and Gmi (Qi). 


however this is likely the result of a significant increase 
in the uncertainty of at early time-slices. Guided 
by the fit-windows at heavier masses we choose to fit 
at later times, however it is certainly possible that we 
simply do not have sufficient statistics and so obtain a 
value for that is suppressed. Another possibility 

is that we are sampling too far into the time-like re¬ 
gion for the VMD hypothesis to hold. Therefore this 
result may suggest that this process is suppressed at large 


TABLE I. Lattice results for the charge-weighted form factor 
Gmi of p^ —>■ 7r^7 for the two available Q^-values for each 
pion mass. 


K 

(GeV) 

g" (GeV^) 

Gmi(Q") 

13700 

0.623 

-0.042(5) 

0.572(12) 



0.078(3) 

0.498(7) 

13727 

0.514 

-0.085(5) 

0.559(13) 



0.070(3) 

0.470(6) 

13754 

0.389 

-0.163(6) 

0.572(15) 



0.066(3) 

0.442(8) 

13770 

0.284 

-0.246(10) 

0.578(35) 



0.054(6) 

0.401(9) 

13781 

0.161 

-0.370(17) 

0.413(64) 



0.050(9) 

0.389(15) 


time-like momentum transfers. As a check of the poten¬ 
tial impact that this may have on the extracted value of 
Gmi( 0), we compare our results with the VMD estimate 
obtained using space-like data and the rho meson mass 
as the monopole mass. This result is shown as the dashed 
line in Fig. [31 In Fig. 0] we show the monopole masses 
extracted from our analysis. We also include the corre¬ 
sponding rho meson masses and find that they compare 
reasonably well at large quark masses, but diverge in the 
light quark regime. Indeed, at the lightest mass A ~ 0. 
This draws the monopole ansatz into question deep in 
the time-like regime. Thus to connect to experiment we 
also work with results in the space-like regime (Q^ > 0) 
and employ VMD. 

We now consider the full hadronic transition form fac¬ 
tor. As was noted earlier, for this transition the quark 
and anti-quark sector contributions are of equal magni¬ 
tude and opposite sign and so the charge weighted form 
factor for p'^ —>■ 7r“''7 is given by 

Gmi{Q^) = IgIi,(Q^) + ^{-Glj.iQ^)) 

Table |T] summarises our lattice results. In Fig. [5] we 
present our results for Gmi (0). We also include the non- 
relativistic quark model (NRQM) expectation and 

the available experimental data. Due to the different 
choice of normalisation for this matrix element, we match 
conventions via the decay width. The relevant expression 
for the decay width using our choice of normalisation Q 
is 

rp^^^ = -a^|GMi(0)P, (10) 

3 nip 

where the photon 3-momentum is evaluated in the rest 
frame for p meson 
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Using this expression we are able to evaluate Gmi( 0) 
for experimental measurements of the decay width. We 
include the PDG average nj as well as the three experi¬ 
mental measurements [37H39j used in its evaluation. For 
the quark model, the choice of normalisation used in Ref. 
[ 3 ^ results in the following expression for the decay width 

As discussed in Ref. [H], by using SU(6) quark and anti¬ 
quark flavour combinations, the sum evaluates to 



Matching with Eq. m, this gives rise to the following 
expression 

Gmi(0) = ^ mpE^ — . 
o e 

For the quark moment, Hud, we use a constituent quark 
mass that varies linearly in 

mud = a + bml, 

with a and b hxed such that the constituent quark has 
a mass of 330 MeV at the physical point and 510 MeV 
at the SU(3) symmetric point, as determined in Ref. [l| 
using the magnetic moments /ip, /i„ and /ta- 

Beginning with our results at the heaviest values of 
m^n-, we find that our lattice data is close to the quark 
model expectation. Furthermore the observed trend in 
the data suggests consistency with increasing quark mass 
where we would expect the quark model expectation to 
hold. The trend in the data also appears consistent with 
the determination of Ref. @ . As we move down to the 
lighter masses there is a clear downwards trend in the lat¬ 
tice data, as was also observed in our previous quenched 
study 0 US- In contrast to this the quark model result 
which shows little variation with varying ttItt. In the light 
quark regime, we do not expect the quark model to nec¬ 
essarily hold for this transition. Unlike heavier systems 
such as bottomiuin and charmoniuin decays, the quarks 
in both the pi and rho meson systems are highly relativis¬ 
tic. Furthermore, the pion itself is a goldstone mode of 
QCD stemming from the underlying chiral symmetry of 
the theory. Th oug h we can in principle treat the quarks 
relativistically [^, , the inability to properly describe 

the chiral behaviour of the pion is a fundamental shortfall 
of all constituent quark models and may certainly lead 
to deviations in the chiral limit. 

Comparison with the experimental determination 
shows a notable deviation, with the lattice data sitting 
around 33% lower than the experimental value. In Fig. [5] 
we also include the values obtained using VMD with the 
space-like data only and find that significant differences 
persist. However, we note that our calculation is incom¬ 
plete. Unlike the elastic form factors for which discon¬ 
nected contributions are necessarily zero [isl . [3l|, such 


FIG. 5. Results for the full hadronic transition moment, 
Gmi(0), extracted using both space and time-like values 
(blue circles) and using the space-like values only and invoking 
the VMD hypothesis (red squares). We include the available 
experimental extractions (orange) ordered from left to right 
by year of publication f37H39| . We also include the PDG aver¬ 
age (red) obtained from this data [I[. The experimental data 
are offset for readability, with the PDG average aligned at 
the physical point indicated by the vertical dashed line. The 
dashed green line is the non-relativistic quark model expec¬ 
tation of [ 3 ^ as discussed herein. We also include the result 
of Shultz et al. (black) which determines this moment for 
a single heavy quark mass. 


contributions are present for meson transitions that in¬ 
volve a change in G-parity [^. This is related to the 
invariance of the QCD action under charge-conjugation. 
For disconnected contributions involving the electromag¬ 
netic current, under charge-conjugation the “G-even” 
two-point function and the “G-odd” disconnected loop 
give rise to a relative sign between the {U} and {U*} 
configurations and so cancel exactly [3l| in the ensemble 
average. For transitions involving a change in G-parity, 
the two-point function is now “G-odd” and so combin¬ 
ing with the disconnected loop gives rise to a common 
sign between the {U} and iC*) configurations result¬ 
ing in a non-zero quantity [9|, ll5|. Furthermore, if one 
neglects disconnected s-quark contributions, the charge 
weight factors between the connected contributions and 
the disconnected contributions are equal 


Gmi 


— n r'con. 

— Hu C^Ml 

2 I 

3 ~ 3 

_ ^con. 

-3 Gmi 



+ (“G'mi ) + qu Gmi + G^n 
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Thus the discrepancy between our results and the ex¬ 
perimental value suggest that disconnected contributions 
are likely to play an important role in fully describing 
this transition. One would also expect such contribu¬ 
tions to become increasing important with decreasing 
quark mass. This expectation complements the obser¬ 
vation that our results are consistent with quark model 










expectations at heavier masses, however deviate as we 
move to light-quark regime. 

Another important aspect that warrants further con¬ 
sideration is the fact that the rho meson is ultimately a 
resonant state. The nature of resonant states on the lat¬ 
tice is signihcantly different to that of continuum due to 
the lack of a continuous distribution of momentum modes 
and so to properly make connection with the continuum 
expectation, one must suitably evolve the lattice deter¬ 
minations to the infinite volume to properly account for 
the differences in the underlying multi-particle interac¬ 
tions. Understanding exactly how to do this in general 
is an area of current interest to the community and only 
very recently has a framework been presented to handle 
1 —>■ 2 body processes required to properly describe this 
transition In fully addressing this aspect, it is 

important to also consider the role of tttt scattering states 
in the p-meson correlation matrix to ensure a complete 
isolation of QCD eigenstates on the finite volume. 

VI. CONCLUDING REMARKS 

In this work we present the first light quark examina¬ 
tion of the radiative decay of the rho meson. The calcu¬ 
lated transition moment, G'mi(O) was found to be con¬ 
sistent with quark model expectations at heavy masses. 
However we have discovered an important quark mass de¬ 
pendence. Our results in the light quark regime sit low in 
comparison with experimental determinations, suggest¬ 


ing important disconnected contributions to this process. 

These results warrant a more comprehensive investiga¬ 
tion of this process. Any future work should aim to fo¬ 
cus on the inclusion of disconnected contributions, multi¬ 
particle contributions and finite volume effects so as to 
allow for the proper evolution of the lattice determina¬ 
tion to the infinite-volume. Understanding the role these 
systematics play in this calculation may provide impor¬ 
tant insights into topical transition amplitudes such as 
N* Nj* and ^ Nj* transitions, and those cen¬ 
tral to the search for exotic hadron states Q. 
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